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The spin-degeneracy of Bloch bands in a crystal can be lifted^ w hen 
spin-orbit (SO) coupling is present and inversion symmetry is absent. In 
two-dimensional electron systems (2DES) spin-degeneracy is lifted by Rashba 
interaction 3 4 terms - symmetry invariants that are scalar products of spin and 
orbital axial vectors. Rashba interactions are symmetry allowed whenever a 
2DES is not invariant under reflections through the plane it occupies. In this 
paper, we use a tight-binding model informed by ab initio electronic structure 
calculations to develop a theory of Rashba splitting in the t 2g bands of the two- 
dimensional electron systems^^ formed at cubic perovskite crystal surfaces and 
interfaces. We find that Rashba splitting in these systems is due to atomic-like 
on-site SO interactions combined with processes in which t 2g electrons change 
orbital character when they hop between metal sites. These processes are absent 
in a cubic environment and are due primarily to polar lattice distortions which 
alter the metal-oxygen-metal bond angle. 

Bulk cubic perovskites have chemical formula ABO3 and the crystal structure illustrated 
in Figure [l^,. The 2DESs in which we are interested are formed from conduction band 
.B-site transition metal <i-orbitals. Because the B site, at the cubic cell center, has octa- 
hedral coordination^ with neighboring oxygen atoms, located at the centers of the cubic 
cell faces, oxygen-metal bonding partially lifts the degeneracy of the d-orbitals, pushing the 
e g = {x 2 — y 2 , 3z 2 — r 2 } orbitals up in energy relative to the t 2g = {yz, zx, xy} orbitals (Fig- 
ure Hp). In the the simplest model of the bulk electronic structure^, the bonding networks 
of the three t2 g bands are decoupled; an o^-orbital on one B site, for example, can hop only 
along the y or x direction through an intermediate p x or p y orbital to an xy-orbital on the 
B site of a neighboring cubic cell. In perovskite 2DESs, the t 2g bands are reconstructed^^ 
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into 2D subbands whose detailed form depends on the bulk band parameters^, the surface 
or interface confinement mechanism, and the dielectric response of the material. A polar 
displacement of A and B atoms relative to the oxygen octahedra occurs in response to the 
confinement electric field; this is the same response that is responsible in some materials (in- 
cluding in particular SrTiOs) for extremely large bulk dielectric constant d 12 * 15 ! At the same 
time atomic-like SO splitting interactions hybridize the three t% g orbitals, which are split in 
the 2DES by differences in their confinement energies. (For (001) 2DESs (assumed below) 
xy orbitals, which have weak bonding along the z-direction, have the lowest confinement 
energy and therefore higher occupancy than {yz, zx} bands.) We explain below how these 
two effects combine to produce a Rashba interaction. 

The Rashba interaction couples an orbital axial vector that is odd under z — > —z to 
spin, and must therefore arise from hopping process that are odd under inversion in the 
x — y plane. We therefore begin by considering a single plane (see Figure ^p) of metal atoms 
and identify the relevant process by using a tight-binding modeP^ for p — d hybridization, 
assigning a hopping amplitude t pd to the process discussed above. Because of the difference 
in parity between p and d orbitals, t p d changes sign when the hopping direction changes 
(Figure [2]). To leading order in t p d virtual hopping via oxygen sites, the Hamiltonian is 
diagonal in the ^-space with eigenenergies: 



e yz = 4t i — 2t i cos (k y a) (1) 
e^x = 4^ - 2^ cos(k x a) (2) 
e xy = 4t 2 — 2t2Cos{k x a) — 2t 2 cos(k y a) (3) 
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Here tip = tp d /A p d where A p d is the splitting between the oxygen p and metal t<i g energy 
levels and the subscripts acknowledge a symmetry allowed difference, ignored below, between 
xy and {yz,zx} hopping amplitudes in the planar environment. Note that the xy band is 
twice as wide as the {yz, zx} bands and lower in energy at the 2D T point. Level repulsion 
from apical oxygens contributes 2t\ to e yZjZX . Because effective metal-to-metal hopping 
amplitudes in this model are independent of hopping direction, they do not produce Rashba 
splitting even when combined with on-site SO terms. 

Rashba interactions are caused by broken mirror symmetry and in particular by the 
associated electric field E perpendicular to the 2DES plane. For t 2e 2DESs, this field both 
polarizes the atomic orbitals and induces a polar lattice displacement. These effects open 
new covalency channels in the metal-oxygen network. In particular there is no hopping in 
the unperturbed system between a metal zx-orbital and an oxygen x-orbital separated along 
the ^/-direction. This is because the x-orbital is even and the zx-orbital odd under reflection 
in the x|/-plane passing through the metal-oxygen bond. When E 7^ 0, the Hamiltonian is no 
longer invariant under this reflection and the hopping process is allowed. If, for example, we 
think about the perturbation as arising from an additional potential —eEz, we can write the 
induced hopping amplitude approximately as Eji, where 71 = (zx, R = 0| — ez\x, R = a/2y) 
(Figure |2^i). At the same time, the electric field will produce forces of opposite sign on metal 
cations and oxygen anions. The induced polarization will change the metal oxygen bond 
angle introducing a non-zero z-component direction cosine n in the bond axis direction. In a 
2-center approximation, this change also gives a non-zero amplitude nt p d for zx to x hopping 
along the y direction. Similar considerations imply an identically induced yz to y hopping 
amplitude along x. (See Figure]^). Including these weak effects, which at leading order 
act only once in the two-step metal-oxygen- metal hopping process, we obtain an additional 
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effective metal-to-metal hopping amplitude that changes sign with hopping direction and 
therefore produces a Rashba effect. The (yz, zx, ^^-representation Rashba Hamiltonian is 



H7 



^ — 2Ur sm(k x a) 

— 2itft sin(k y a) 

y 2Ur sm(k x a) 2it R sm(k y a) J 



(4) 



where the Rashba interaction strength parameter tn = (jit pd E) / A pd + (nti d ) / /S. pd . When 
combined with the an atomic-like bulk SO interaction 1 12 1 14 1 , described in the Supplementary 
Information, H^ 9 leads to Rashba splitting in the t^g bands. We remark that broken mirror 
plane symmetry also introduces other covalent bonding channels, but these do not contribute 
to the Rashba effect. We note that, a surface metal atom in a BO2 terminated perovskite 
is not octahedrally coordinated. This absence of local inversion symmetry and the decrease 
in level repulsion with neighboring oxygen atoms mixes e g and ti g orbitals at the surface. 
When this mixing is strong a more elaborate theory of Rashba SO coupling is required. 

In general t<i g 2DESs will be spread over many coupled metal layers, and the Rashba 
Hamiltonian H 1 ^ will act within each layer with a layer-dependent coupling constant t R . For 
the extreme case of a single-layer t2 g 2DES, the xy- band will be pulled below the {yz, zx}- 
bands by differential confinement effects. In this case we can derive a simple effective Rashba 
Hamiltonian which acts within the xy subspace. To do so, we define 5 as the energy scale 
which splits the xy and {yz, zx} bands at the T point. Allowing virtual transitions to the 
{yz, zx} manifold due to orbital/lattice polarization (He), and bulk SO effects (Hso), we 
find the part of the Hamiltonian linear in electric field is given at small k by, 

H x R y = e xy (k)-aa-(kxz) (5) 
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where a = AAsotRa/(35). 

To support our theory of the Rashba effect we have carried out an ab initio study of 
a typical t 2g 2DES. To simplify the comparison we examined the case of a single (001) 
BO2 plane and studied the influences of z-direction external electric fields and oxygen- 
metal sublattice relative displacements separately. Because we expect Rashba splitting to 
be proportional to A$o we use the 5d transition metal Halfnium (Hf) as the B atom. To 
minimize the mixing between t 2g and e g bands apical oxygen atoms have been included in 
our study to maintain octahedral coordination of the metal sites and maximize crystal field 
splitting. Figure [3^i shows the band structure of a Hf perovskite plane with ideal atomic 
positions in the absence of an applied external electric field when spin-orbit interactions are 
neglected. At the zone center, the xy band has a lower energy than the {yz, zx} bands 
as expected in t 2g 2DESs. The strength of the Rashba hopping processes can be read off 
the band structure by identifying the avoided crossing which occurs between an xz or yz 
band and the small mass xy band along a large mass direction in momentum space when an 
electric field is present. We find that even for an extremely large electric field, Q.leV/K , the 
level repulsion (Eqn. [4]) at the crossing is very small. Figure |3}d shows the band structure 
changes when SO coupling is included. Note the expected confinement-induced t 2g manifold 
degeneracy lifting at the T point. On the scale of this figure the Rashba splitting is too 
small to be visible. Figure [3^ plots the spin splitting, which is largest near the band crossing 
and reaches a maximum value of ~ 5meV, as a function of k, and compares it with the 
splitting predicted by our theory when the value of tf> is fit to the ab initio bands calculated 
in the absence of spin-orbit coupling (see Supplementary Table I for tight-binding fitting 
parameters). 

Figure [4] reports the corresponding results obtained for the case of a polar lattice displace- 
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ment. Figure |4^l shows the bandstructure for a Hf perovskite plane with a displacement of 
0.2 A. The Hf atom displacement was chosen to emulate the interface atomic configuration. 
According to previous first principles calculation in LaTiOa/SrTiOs interface^, Ti atoms 
at the interface are displaced out of the Ti02 plane by « 4% of the lattice constant. For 
illustrative purposes we report results for a similar Hf atom displacement of ~ 5% of a lattice 
constant. 

The yz, xz-xy avoided crossing is now easily visible. After SO coupling is included, a clear 
Rashba splitting is visible (Figure |4|d,c) that is an order of magnitude larger than for the 
orbital polarization case. Note that our t 2 g only model underestimates the zx band splitting 
in both cases, particularly close to the band center. We ascribe this to the proximity of the 
lower e g level that is visible in the band plots and neglected in our theory. 

Our theory of Rashba interactions in t2 9 2DESs is consistent with experimental evidence 
for strong spin-orbit interactions at interfaces between polar and non-polar perovskites^^Q], 
for example the SrTiOa/LaAlOs interface, and in surface 2DESs induced by the very strong 
electric fields applied by ionic liquid gates. It also suggests that Rashba interactions will tend 
to be stronger in materials which are more easily polarized. In this sense SrTiOs has a poten- 
tial for relatively strong spin-orbit interactions even though Ti is a 3D material. t% g 2DESs 
with strong lattice polarizabilities may^l contain both weakly-confined orbitals responsible 
for superconductivity and strongly-confined orbitals responsibility for magnetismP^. If so, 
our theory implies that magneto-transport properties in these materials will be strongly 
sensitive to local lattice polarization at the surface or interface. We believe that this work 
provides a starting point for the interpretation of heretofore unexplained magneto-transport 
phenomen a! 22 * 23 ! 

Methods: The first principles calculations were based on the density functional the- 
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ory and carried out using the Vienna Ab Initio Simulation Package^. We used projector- 
augmented wave pseudopotentials and the generalized gradient approximation exchange- 
correlation functional of Perdew, Burke and ErnzerhoP^. The supercell contained a three 
atom HfC>2 layer with two oxygen atoms located directly above and below the Hf atom. The 
molecular layers were separated by 20 A vacuum. The plane-wave energy cutoff was set to 
500 eV. We employed a 8 x 8 x 2 /c-point sampling to achieve electronic convergence. 
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FIG. 1: Perovskite crystal structure, a, Bulk cubic unit cell with the A atom in blue, B atom 
in red, and the oxygen in white, b, Splitting of atomic d-orbitals into e g and tig manifolds, c, 
Single BO2 plane with one unit cell, with area a 2 , in the boxed region. 
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FIG. 2: Bonding network along the y-axis with an electric field, a, Orbital polarization. 
Bonding between zx and xy on neighboring metal atoms through p x orbitals. The positive and 
negative lobes of the orbital functions are represented in blue and red, respectively, b, Lattice Po- 
larization. Displacement of the metal (light blue plane) and oxygen (light purple plane) sublattices 
in an electric field. 



12 



TT/4 ft/2 3TT./4 TT 

k x [l/a] 



FIG. 3: Orbital polarization changes to t2 g bandstructure. a, t<i g bandstructure in the 
absence of an electric field and SO coupling. The results of our model are shown as solid lines 
while the simulation is shown as dotted lines. The band order (at the T point) is xy (blue), yz 
(red), zx (green), and x 2 — y 2 (purple). The inset shows the onset of an avoided crossing in the 
presence of an electric field, b, tig band structure with SO coupling and orbital polarization, c, 
Comparison of the orbital polarization part of the Rashba splitting in the tig space. 
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FIG. 4: Lattice polarization changes to t2 g bandstructure. a, t2 g bandstructure in with 
lattice displacement and no SO coupling. The results of our model are shown as solid lines while 
the simulation is shown as dotted lines. The band order (at the T point) is xy (blue), yz (red), zx 
(green), and x 2 — y 2 (purple), b, ti g band structure with SO coupling and lattice displacement, c, 
Comparison of the lattice mediated Rashba splitting in the t2 g space. 



14 



Supplementary Information 

The explicit form of the atomic-like SO coupling Hamiltonian^^ mentioned in the main 
text is given in the (yz f, zx f, xy f, yz I, zx I, xy \) representation by, 
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Breaking the mirror symmetry of the 2DES allows additional covalent bonding terms - not 
present in the unperturbed cubic system. Those described in the main text are associated 
with the Rashba effect, but others new processes are allowed. These are summarized in the 
Table II. In the unpolarized system, the p x orbital is not involved in the bonding along the 
x direction. For this reason 72 processes, which couples zx and p x along the x direction, can 
only contribute to second order in E. The corresponding statement is also true for lattice 
displacements. 

73 increases or decreases the bonding strength with the apical oxygens. If both apical 
oxygen atoms are present, it can only contribute a term in the Hamiltonian that is quadratic 
in electric field. At the surface of a BO2 terminated perovskite, there is a missing apical 
oxygen. Therefor, within the ^ space and to first order in the electric field, the following 
term must be added to the Hamiltonian of the surface layer: 



' A s 0^ 



V 



/ 



In the above equation, A< 



Hs A s • (S2) 


E^tpdn I A p d and may be positive or negative. Because 
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Tight-Binding Model Parameters 


Lattice Constant 
SO splitting 


a 


4.05 A 
0.340 eV 


Orbital Polarization 


t 2 

tR 

a 


0.41 eV 
0.51 eV 
0.0014 eV 
0.02 eV 


Lattice Polarization 


ti 
t 2 
t R 
a 


0.30 eV 
0.48 eV 
0.045 eV 
0.02 eV 



TABLE I: Tight-binding parameters used in fitting ab initio results for BfC>2 plane. 



Tight-Binding 


I Matrix Elements 


Lattice Polarization Effects 


(xy, R = 0\AU\{x, y, z), R = ±%x + n%z) 
(xy, R = 0\AU\{x, y, z), R = ±\y + np) 
(yz, R = 0| AU\{x, y, z},R = ±%x + np) 
(yz, R = 0\AU\{x, y, z},R = ±%y + n%z) 


{0, ±t pdw , 0} 
{±t pdn , 0, 0} 
{0,nt pdn ,0} 
{Q,nt pdw ,0} 


Orbital Polarization Effects 


(xy,R = 
(xy, R = 
(yz, R = 
(yz,R = 


— ez\{x, y,z},R = ±|x) 

- ez\{x,y,z},R = ±|y) 

— ez\{x, y, z}, R = ±|x) 

- ez\{x,y,z},R = ±fy) 


{0,0,0} 
{0,0,0} 
{0,7i,0} 
{0,72,0} 


Apical Oxygen Atoms 


(xy, R = 
(yz,R = 


— * 

- ez\{x,y,z},R = ±p) 

- ez\{x,y,z},R = ±p) 


{0,0,0} 
{0, ±73,0} 



TABLE II: Tight-binding matrix elements for metal tig and oxygen p-orbitals. The matrix elements 
for the zx orbitals can be derived from the yz entries, by symmetry. 

there is no periodicity in the z- direction, this term is independent of k. In addition, the 
lack of octahedral coordination at the surface can lead to significant t 2 g/e g hybridization. 

In addition to the bonding changes to the cubic system mentioned above, if the structural 
distortions are more complicated (e.g. there is an in-plane twisting of the oxygen octahedra) 
some a bonding between the oxygen p and metal t 2g orbitals will also be allowed. Provided 
the polar displacement is still present, these a bonding channels can also contribute to the 
Rashba effect. 
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